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ABSTRACT 


Since a nonlinear feedback control system may possess more than one 
type of forced oscillations, it is highly desirable to investigate the 
type of forced oscillations which can occur when the nonlinear restoring 
force function is of a specific form, 

This paper is to present a harmonic linearization method for finding 
the existence of forced oscillations and response curve characteristics 
of a nonlinear feedback control system by means of finding the restoring 
force function of the nonlinearity and using the harmonic balance and an 
iteration method for investigating the conditions for one type of forced 
oscillations exhibited. 

The existence conditions for fundamental frequency, sub-harmonic of 
2nd order and 3rd order forced oscillations of a second order feedback 
control system are investigated; also the fundamental frequency forced 
oscillation for a higher order system and the jump resonance frequencies 
of response curve are investigated, and a general expression for the equi- 
valent gain of the nonlinearity has been developed. 

The author wishes to express his appreciation for the assistance and 
encouragement given hy Dr. George Julius Thaler of the U. S. Naval Post- 


graduate School in this investigation. 
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CHAPTER I 
INT RODUCT ION 

One of the major branches of servo mechanism analysis for nonlinear 
control systems is the study and methods of prevention of continuous oscil- 
lation. For a linear control system, the output may oscillate with con- 
tinuous or increasing amplitude with no input signal applied, in which case, 
the system is said to be unstable, or they may have oscillations which die 
away. In both cases this transient response is the characteristic of the 
system itself and the conditions for divergent, continuous or decaying 
oscillation do not depend on the form or magnitude of the input signal. 

When a nonlinear element is present, this independence of the input 
signal ro longer holds. The characteristic performance of the nonlinear 
control system will depend on both, i.e., forcing function and the 
characteristic of the nonlinear element. 

A nonlinear control system, under suitable conditions may exhibit 
steady oscillations in which the main component has a frequency which is 
dependent on the frequency of the forcing function. This type of oscilla- 
tion is called "Forced Oscillation". A forced oscillation for which the 
frequency is a fraction of the forcing function frequency is called "Sub- 
harmonic Forced Oscillation’. A forced oscillation for which the frequency 
is multiple of the forcing frequency is called "Super Harmonic Forced Oscil- 
lation", 

Since a nonlinear control system may possess more than one type of 
forced oscillations, it is highly desirable to investigate the type of 
forced oscillations which can occur when the nonlinear restoring force 
function is of a specific form. 


Some methods of investigation for nonlinear forced oscillation have 





been discussed by J. C. west! in his book "Analytical Techniques for 
Nonlinear Control System (1960). A graphic method for investigating 
forced oscillation of ''ON OFF" control system has been developed by Hamel 
(French) and Tsypkin (ussR)?. Some work for subharmonic oscillations of 
a typical restoring force function has been done by C. A. Ludeke and 
William Paae in their paper 1959, and also by Ogata in his PhD pieeis *, 
June 1956. 

There are some other authors for investigation of the vibration of a 
mechanical system, such as N, Wimersky”. J. J. sueuer, and Y. H. ar in 
their books of ‘Introduction to Nonlinear Mechanics";''Nonlinear Vibrations" 
and "Analysis and Control of Nonlinear System'' respectively. 

The purpose of this paper is to present a harmonic linearization 
fiechod for finding the existence of forced oscillations and response 
curve characteristics of a nonlinear feedback control system, by means of 
finding the restoring force function of the Pomiimemtity’? so and using an 
iteration method for investigating the condition for one type of forced 
oscillations exhibited. A response curve and "Jump Resonance" can be 
also investigated, 

The emphasis of this paper is placed on showing the general approach 
of the method, therefore, second order nonlinear feedback control systems 
are used as examples. However, a higher order system will be discussed in 
this paper also, 

By the harmonic linearization and iteration method, a general equa- 
tion for the equivalent gain of an odd function nonlinear element will be 
veloped: ', It is in terms of the amplitude and the frequency of input. 

There are six charters in the main body of this paper. The first 


chapter and the last chapter are a general introduction and conclusion 





respectively. Chapter two is a general description of oscillation of 
feedback control systems and an investigation of the restoring force func- 
tion of the nonlinear element. Chapter three investigates the existence of 
fundamental forced oscillations and the jump phenomena in the nonlinear 
feedback control system. Chapter four is a sub harmonic oscillation in- 
vestigation. Chapter five investigates the existence of forced oscillation 


of a higher order feedback control system, 





CHAPTER II 

GENERAL DISCUSSION OF OSCILLATION IN FEEDBACK CONTROL SYSTEMS 
2-1 General Description: 

For either a linear or nonlinear feedback control system, most re- 
quirements for control system design specify a minimum time response, high 
accuracy, and high stability. It is clearly in order to have a minimum 
time response and high accuracy for a feedback control system, but first 
of all, the system should be stable. 

In general, two types of oscillations may be exhibited in a feedback 
control system, one is free oscillation, i.e., an autonomous system, This 
is a system which is left by itself with its outside source of energy sudden- 
ly cut off or abruptly changed in amplitude, In such a system, free oscilla- 
tions occur as the system tends to equilibrium after reaching a state of un- 
equilibrium. The oscillation frequency of this type depends only on the 
characteristic of system itself. Ordinarily, a free oscillation is a damp- 
ed oscillation; as the time increases the amplitude decreases and damps out 
in steady state. For an ideal case, if it is a system without damping, i.e., 

Xt = 0, the oscillation will be continuous. 

The other type of oscillation is a forced oscillation, i.e. a nonauto- 
nomous system. This is a system that is acted upon by an outside source of 
energy, a force function. The force function may be a constant; a periodic 
function or any other function of time. 

The analysis and methods for preventing oscillations of a feedback 
control system are concerned directly or indirectly with the characteristic 
of system, i.e., the differential equations of the system. If the system 
is linear a frequency response of the system can be directly solved from 
the differential equation, furthermore, the principle of superposition af 


be applied too. On the other hand, if the system is nonlinear, either with 





a nonlinear damping or nonlinear restoring force function, then the basic 
tools for the analysis of a linear system are no longer valid. 
2-2 Oscillation in Linear Systems: 

There are two types of methods for analysis and design of linear feed- 
back control systems; one is graphic method, the other is an analytical 
method. The best known of graphic methods consist of Bode diagram; Polar 
plot; Nichols diagram; and Root locus. If we know the transfer function of 
the system, any of the above methods can be used for solving the problem of 
stability. The analytical method is a mathematical analysis for solving the 
differential equation, and plot the curves of the solved equation. 

First consider a simple linear feedback control system, as shown in 


Fig. 2-1. The forward element has a transfer function G(s), and a feedback 


Fig. 2-1 Block Diagram of a Linear System 
element transfer function H(s). The feedback signal is H(s)® (s) and it 
is assumed that this is subtracted from the input signal 6, (s) to form 
the control signal e(s); thus: 
e(s) = 9,(s) - H(s) 2 (s) (2-1) 
The relation between the control signal and reference signal in the 
Laplace transformation form is: 


G(s 


a) 
&. (s) - 1 + H(s)G(s) &,(s) (2-2) 





The solution of equation (2-2) for a particular input can be obtained 
by Laplace and Heaviside methods, This involves the determination of the 
roots of equation: 

1 + H(s)G(s) = 0 (2-3) 
In order to split the transfer function into the sum of partial fractions 
of standard form for which time solutions are known, If equation (2-3) is 
ann th order polynominal in s, there will be n roots say S, to s? such 


1 


that (s- s) is a factor of the polynominal. Special cases arise in which 


some of roots are identical, and if m roots of SH occurs, (s - $,)" is a 
factor of the polynominal. 
In general, the roots 7 will be a complex number, as: 
S. mz ev a ct (2-4) 


in which of and a are both real numbers, and since the coefficients of 
the polynominal of equation (2-3) are all necessarily real, then the com- 
plex roots should be in conjugate pairs, Thus: 
Ce ees ae CS Ke) 
will be a factor of the polynominal, 
In determining the transient response to a forcing function input, 
these various factors become the denominator of the partial fraction re- 
presentation of the transfer function equation (2-2), The term of KR /(s -S_) 
gives a time solution containing expoential term: 
AG oe hace Dt 
The complete solution is the sum of all such terms, 
The complex conjugate pairs give rise to an oscillatory term: 
e“r® sintw t+ ¢ ) 
If the system has convergent response, all of these terms in the time solu- 
tion must converge and hence all roots S, to sf must have negative real 


l 


parts. If any of these roots has a zero real part, then a continuous 


6 





oscillation is produced. If any of the complex roots has a positive real 
part, then divergent oscillation will be the result. 

To determine the stability of a system by using the graphic method, it 
is not necessary to find the roots of the characteristic equation. If the 
equation is known algebrically, the Routh method can be applied and also 
by the Nyquist and Bode diagram method a graphical solution can be obtained. 

Nyquist criteron graphical analysis is the most useful method to solve 
the stability problem of control systems. Rewriting equation (2-3) in the 
form of: . 

G(s)H(s) =-1 (2-5) 
The roots of equation of (2-3) are particular values of s which satisfy 
equation (2-5). Thus if the equation of G(s)H(s) can be mapped for continu- 
ously variable values of s, then those values for which G(s)H(s) = -1 can 
be determined. 
Assume a particular roots of S apt ON: ju can be represented by a 


particular value of G(s)H(s), i.e. 


G(< + Ju HCL | + jiu » = R(/,, We) + jl ,e! ) (2-6) 
This equation represents a vector in I vs R plane, the magnitude of vector: 
We! lw (2-7) 

and the direction by an angle: 
¢ = Tan 'T/R (2-8) 


For a constant < _, the value of M will be a locus in I-R plane as 
(2 varies. This process can be repeated for different values of :4 r_ to 
obtain a family of curves as shown in Fig. 2-2. 

This steady state frequency response locus of the open loop system 
can be obtained experimentally from the system by measuring gain and phase 
for a sinusodial input of varying frequency. It can be shown that is the 
Steady state frequency response locus passes to the right side of (<1, 0) 


ri 





Fig. 2-2 Locus of an Open Loop Transfer 
Function G(s)H(s) 


point, this system will be stable, on the other hand, if it passes 
through or to the left side of (-1, 0) point, the system will be unstable. 
This criterion is due to Nyquist. A typical curve is shown in Fig. 2-3. 

In the extension to nonlinear control systems, or more complex multi- 
loop systems, the single loop criterion is sufficient for the majority of 
needs where G(s) can represent the over all transfer function of several 
loops provided that G(s) is itself stable. 

zs 
— 


ae 


c) | R 


(=6, 
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/ 
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Fig. 2-3 Nyquist Criterion for Three 
Characteristic States 





2-3 Forced Oscillations in a Linear System: 
Assume a second order linear system block diagram is shown in Fig. 


2-1; in which: 
K 


G(s) = te C8) (2-9) 
and: 
H(s) # 1 (2-10) 
From equation (2-2), the differential equation of system becomes: 
@ +76 + KO = KO (2-11) 
c Cc Cc < 
Let the forcing function: 
@ (t) «x F Coswt (2-12) 
Equation (2-11) becomes: 
9 vt LO Ut CKO = 6C Coswt (2-13) 
Cc Cc Cc 


in which C is a constant of value of FK. 

The solution of equation (2-13) consists of the sum of the solution 
of the homogeneous equation, (i.e. the free oscillation of the system) and 
the solution of the non-homogeneous equation. Let the solution of homo- 
geneous equation: 


Q(t) = ae on Cosw fist + C, Sin Wns 17 t) (2-14) 


in which: 
ne AK 
eo me Jy 
f # Cos\ 


The complete solution of equation (2-13) will be: 


Sinw myc) 0 ~15) 


i oe Wo jer Ai 


@(t) ze FG Cosi Jjl-f t+, 


Equation (2-15) is obtained by a superposition of the free oscillation 
and the forced oscillation which varies from the action of external force: 


The frequency of forced oscillation is the same as that of the external 





force, the magnitude of the forced oscillation is given by: 


C 
: z = ga* ay 
{ Ff “71 ° by ¥, t ef dx ? Le 


The value of phase shift @ related to the external force is given 


by: 


Cos® = a a ee : Oey) 


. yx os ‘ 
, (uw 4 tu *) + Tf wo AW 





— we —_ (2-18) 


— ty ee ¢ fb PEP a 


~~) 





Sin@ =. 


; 
As 
 ) 


I 


In the case of positive denpine . i.e. <> 0, it is clear from equation 
(2-15) that after a sufficiently long period of time, the free oscillation 
is damped out and only the forced oscillation would be observed. 

In the case of no damping, i.e. «< = 0, the phase shift 0 is seen 
from equation (2-17) and (2-18) to be zero for W<4» and Zz for w>u,. In 
other words, the forced oscillation is in phase with the external force if 
the forced oscillation frequency is less than the free oscillation fre- 
quency, and is 180 degrees out of phase with the external force when ~’ is 
greater than w.. 


In the case of .< = 0, equation (2-15) becomes: 


= Y i 
Q(t) (C, Cosi) t + C, Sin i) t) 


+ C Cos(iit + @) (2219) 


lua — w| 
There will be a superposition oscillation of two frequencies, if the value 
of « +w~ one of which is the natural frequency, and the other is the 
frequency of external force. In the case of w+, the free and the forced 
oscillations have the same frequency, and the solution of equation (2-13) 
will be found: 


C 
6 (t) =C, Coswt +C, Sinwt - zt Sinwt (2-20) 


1 2 2 
In this case, the component due to the external force is no longer periodic, 
it is oscillatory with an amplitude that increases linearly with time. At 
this condition, it is a resonance phenomenon. It is very important in de- 


sign to avoid this phenomenon in a system with a periodic force function, 
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In the case of a damped system, i.e., « » 0, from equation (2-16) for 
the steady state, the amplitude of forced oscillation is always finite. It 
is possible to give, by the use of dimensionless variable, 4 more general 
significance to the expression for the amplitude of forced oscillation, let 


the amplitude: 


C 
|H| = MC me (2-21) 
Therefore equation (2-16) becomes: 
Mom ao (2-22) 


(0-77 Peep)" ( ce) 
Where M is defined to be the magnification factor. The extreme values for 


2 ja 
M are attained for “= 0 and ( “/ ck, x ~#l-2y. %f1l-2/fP=<0, 


there is a maximum for «#0; if 1- 2,/ - 0, and«>Q, there is a maxi- 
mum for w7%'. = 1, -2P° and a minimum for ~)= 0. For small values of damp- 
ing coefficient, the frequency which produced the maximum amplitude is very 
nearly the natural frequency of the system, Fig. 2-4 shows the forced oscil- 
lation response curves of a linear system as a function of /W. with 
various values of / . 
2-4 Oscillations in Nonlinear Feedback Control System: 

As mentioned in section 2-1 of this chapter, when a nonlinearity is 
present in a feedback control system, the properties of proportionality 
and superposition are not valid. In the linear system, these properties 
involve the existence of a transfer function and of characteristic fre- 
quencies pfoper to the system (frequencies at which the system tends to 
oscillate with an amplitude depending on initial conditions); in the case 
of nonlinear systems, the amplitude, like frequency, can depend at one and 
the same time on both the initial conditions and the system itself. Some- 
times, as in the case of limit cycle both frequency and amplitude of input 


are characteristic of the system and independent of of initial conditions. 
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In many cases the analysis and design of such systems using linear 
theory may produce an excellent result. This general procedure is a 
type of linearization procedure, i.e., the actual nonlinear system is re- 
placed by a linear system, which approximates it, and the analysis and de- 
sign techniques are applied to the linear equivalent. The most useful method 
for nonlinear linearization is that of the 'Describing Function" method, It 
is actually a first harmonic approximation method, 

The method for linearization used by this paper is a multi-harmonic 
linearization and iteration method; first assume a harmonic solution for 
the , differential equation and insert it in the differential equation, 
and compare the coefficients of the same order harmonic terms, A detailed 
procedure will be discussed in the next chapter. 

As previously explained, a nonlinear feedback control system may posses 
some type of oscillation. Which type will exist depends on the restoring 
force function of nonlinear element and the amplitude of forcing function. 
Hence, before investigating the existence of forced oscillations, the re- 
storing forcing function of the nonlinearity should be investigated. 

2-5 Characteristic of Nonlinear Elements: 

Nonlinear elements which may be seen in feedback control systems can be 
represented by a combination of four fundamental concepts of nonlinearity, 
which are: (1) Saturation; (2) Variable gain; (3) Dead zone; (4) Hysteresis. 
These basic notions of nonlinearity can, in combination with one another, 
result in almost all types of nonlinearity in practice, as shown in Fig. 2-5. 
For example an ideal relay nonlinearity can be seen as a saturation nonline- 
arity which linear part with a value of infinite slope. 

The operational characteristic of a nonlinear element which operates in 
a control system depends on the operating condition, just as a vacuum tube 


for which the operating characteristic depends on the operating point of its 
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(b) Variable Gain 





fc) Hysteresis 


or Nonlincar 





characteristic curves, As a saturation nonlinearity element with a small 
forcing function input, it operates linearly. On the other hand, if the 
forcing function is very large and with a high frequency, the actual opera- 
ting characteristic {is like an ideal relay. Hence, the actual operating 
characteristic of a nonlinear element may be said to be a function of the 
amplitude and frequency of the forcing function. 
2-6 Restoring Force Function Investigation of a Nonlinear Element: 
Consider the block diagram of Fig. 2-6, in which the output is a func- 


tion of input, usually this function in a control system is called the 





| Nonlinear 


Element 


ee De 
e] 


ere 8 pee ee 


Fig. 2-6 Basic Relation Between Input 
and Output of a Nonlinearty 


"Restoring Force Function". 
Assume the input is: 
ea E(t) (2-23) 
The output will be a function of input, as: 
e, = £(E) (2-24) 
From the numerical analysis and curve ieee process, the function 
of £(E) can be expressed: 


£(E) = ata E+a Ee +a B74 WSs se sibs et aE 


1 2 3 
l/n 


ier + bE!” + cocceeeceseevve00t bE 


+ bE +b (252) 
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Equation (2-25) can be in the general form: 


a Ei + <=: bE (2-25a) 


ECE) % oto . 
a | i'm 


for this paper, the nonlinear characteristic is considered symmetrical to 
the original point, that means: 

£(-E) = -f(E) (2-26) 
It is an odd function, the coefficients of even order of equation (2-25) are 
zero, it becomes: 


f(E) = ao + 2 E + a gE? + a E> Ht ccecnceeececcsecves 


1 3 5 
+ bE + bE! ae bee! Ft cccsccccccceeseseseoec (2-27) 


For the saturation and variable gain nonlinearity, J. C. west! uses 


the restoring force function as a form: 


£(E) = aE + aE” (2-28) 


in which ay is greater or equal to zero, and a, 


or equal to zero corresponding to the case of a "Hard Spring", "Soft Spring", 


is either greater, or less, 


and "Linear Spring" saturation nonlinearity characteristic respectively. 


A set curves for all values of a. is shown in Fig. 2-7. 


JE? as 
| Gy, O 
/ vA 
: Lo &0 
va 


Fig. 2-7 Characteristic curves of Eq. (2-28) 
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In order to generalize for the condition of saturation nonlinearity 
with dead zone, equation (2-28) can be written with a constant as: 


f(E) = -a Sign E+a,E +a E> (2-29) 
oO 1 3 
if the dead zone is less than one, the characteristic curve can be shown 


in Fig. 2-8. 


Fig. 2-8 Characteristic Curves of Eq. (2-29). 

As discussed previously, when a saturation nonlinearity is operating 
with a larger amplitude than the saturation voltage and a high frequency, 
the operating characteristic is like an ideal relay. On the other hand, 
under some suitable conditions a relay nonlinearity control system may ex- 
hibit some sub-harmonic forced oscillations, as proved by the paper of A. M. 
Hopkin and K. Ozate +, This result is the same as a saturation nonlinearity. 

Actually, the operation of a relay is not a perfect discontinuous 
characteristic element, when the relay is operating from static to closed, 
there will be a little amount of time delay from starting to close cone 
tacts. Therefore, the operating characteristic of an ideal relay may be 
considered a saturation nonlinearity with a very large slope of the linear 


part, a sketch of it is shown in Fig. 2-9. 
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Fig. 2-9 Operating Characteristic of 
an Ideal Relay 


From the above discussion, and for a proper value of constants, the 


restoring forcing function can be expressed as: 


3 5 
£(E) = aie - a8 + aE (2=30) 


Similarly, if we consider the dead zone the equation of (2-30) can 
be expressed: 


. ) 5 
£(E) = "a, SignE + aE + aE 7 ack (2-31) 


M. J. Mees investigates the restoring forcing function of an ideal 


relay nonlinearity by a fifth root or cube root of the input as; 


E(E) = bee! (2-32) 


Or; 


E(E) = be! (2-32) 


A more detailed discussion about frequency response by using the re- 


Storing forcing function will be covered in the next two chapters. 
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CHAPTER IIT 
FUNDAMENTAL FREQUENCY FORCED OSCILLATION IN NONLINEAR CONTROL SYSTEMS 
3-1 General Description: 

Consider a system shown in Fig. 3-1, in which G,(s) is the transfer 
function of a controller element. G, (s) is the transfer function of a 
motor and gear element, G, is a nonlinear element. H(s) denotes a transfer 
function of a feedback network. 

When a periodic forcing function is applied to the input, the output 
of the system may be or may not be a periodic function. [f it is a peri- 
odic output and the frequency is the same as the frequency of the input 
forcing function, it is said that the system is operating in forced oscil- 


lation with the fundamental frequency. 





aaa: | 
| | | | 
————— j1¢>) an 4 


ee 


Fig. 3-1 Block Diagram of a Feedback 
Control System 


There are some methods to represent the forced oscillation in a non- 
linear system, If we keep an input forcing function with constant frequency 
and vary the amplitude of input, a response curve can be shown in the "In- 
put vs Error" plane, Fig. 3-2, in which there is a cut off amplitude; that 
means when the input is smaller than ees there is no forced oscillation. 


The other representation is with the input amplitude constant vary 
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Fig. 3-2 Forced Oscillation Represented 
in the E vs 0. Plane 


the frequency of input. It is shown in "Error vs Frequency" plane. Fig. 


3-2a is an example, it is actually a closed loop frequency response of the 


system, 


Mh 





Fig. 3-2a Forced Oscillation Represented 
in the E vs 0 Plane 


Either represented in E vs 0. or E vs c’ plane, both of them can be 


represented in a phase plane (E vs E). A typical phase plane of forced 


oscillation from the computer is shown in Fig. 3-3, in which it is represent- 
ed by an ellipse. The size of it is changed with the frequency. 
3-2 Basic Equation for Forced Oscillation of 2nd Order Nonlinear Feedback 


Control System: 
Consider a 2nd order feedback control system shown in Fig. 3-4: 
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Fig. 3-4 Block Diagram of a 2nd Order 
Feedback Control System [ 


The output equation in Laplace form: 


e(e) « KE (301) 


s(S+.4 ) | 


in the differential form: 


8. + a 8. = Kf (E) (3-2) 
Where: 

@(s) = 8 (s) - E(s) (3-3) 
Equation (3-2) becomes: 

E +4E + KE£(E) = 8 + 6. (3-4) 
Assume the input forcing function: 

Q(t) xe F Cos(wt + 8) (3-5) 


in which @ is the phase difference between the error signal and forcing 


function of input. Inserting equation (3-5) in to equation (3-4): 


BE +E + KE(E) = -ACos(wt + 0) —BSin(wt +0) (3-6) 


or: 
E +26 + KE(E) = C Cos(ut+x+4+0 +4 9) (3-6a) 
Where : 
A= WF (3-7) 
Bazxtu F 7 (3-8) 
[ 2 2 -,\ |e 
C= JA + 8B BRB FWWK +? (3-9) 
G = wear 'B/A = ian wae (3-10) 
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Equation (3-6) or (3-6a) is a basic differential equation of a 2nd order 
nonlinear feedback control system, in which f(E) is the nonlinear restor- 
ing force function, and K is a constant of system. 
3-3 Forced Oscillation; Saturation Nonlinearity Without Damping: 
Consider a system without damping, i.e.-< = 0, equation (3-6) becomes: 
E + KE(E) = -A Cos(wt + 6) (3211) 


Let the restoring force function: 


f(E) = aE + aE (2-28) 
Equation (3-11): 
E * K(a,E + aE) = -A Cos(ut + @) (3-12) 
750 


By the harmonic iteration method , the solution of equation (3-12) can 


be expressed as a series of odd order periodic harmonics: 


e(t) = Ey Coswt + E, Cos3 vt + Eg COStL © sees (3-13) 


y? Ey) ES ooceeee are the amplitude of harmonics of e(t). 


If we consider only the fundamental frequency forced oscillation, we 


in which E E 


can assume one solution is: 


e(t) = Eq Coswt (3-14) 


Where Ey is the fundamental frequency amplitude of error signal, w is the 


frequency of forcing function and forced oscillation; is to be determined, 


Since: 
E =- WE, Cost (3-15) 
E s- i Ey Sinwt (3-16) 
3 1 3 2.4 4 
E = rn Ey (3 CosWt + Cos3Ut) (3-17) 
Inserting equation (3-14) into (3-17) in equation (3-12): 
(- we + Ka,E, + 2 ka E a Cos wt re aE scos3 Wt (3-18) 
1 | Ma sas 80) | 4 3°71 


= -A (Cosivt Cos@ + Sina/t Siné@) 
neglecting the higher order harmonic terms and equating the coefficients 
of CosWt and Sint, thus: 


23 





(10, + Ka E, +2 Ka,E1>) = = ACos@ (3-19) 


ASinwt Sin@ = 0 (3-20) 
Equation (3-20), the term of ASin vt is not always equal to zero, the term 
of Sin6 should be zero. Hence the value of @ will ke either zero or nx, 


Inserting equation (3-7) into equation (3-19): 


Ka.E,+3 Ka.E,3 
|e 8 1 
, 4 1/2 Seva 
Ey —- F 
For @ = 0. and: 3 3 ps 
— tf? 
er (aie ea! laid) (3222) 
_ Ey + F 


for @ = x, 
As a check, in the case of F = 0, i.e. no forcing function input, equa- 


tions (3-21) and (3-22) give the free oscillation frequency of system: 


—- 3 2% 
‘ee i — } = 
we JK (a, + 7 aE) (3-23) 
For the linear case, i.e., a, = 0, the natural frequency 
1/7 
wn = (a,) /- (3-24) 


It is exactly the same as we discussed in the linear system, 
Recall equations (3-21) and (3-22); the value of uw should be real, 
that means: 
ay +t aE, > 0 (3-25) 
for both the phase angle is either zero or x. And : 
Ey vot (3-26) 
for the phase angle is zero only. 
Equations (3-25) and (3-26) are the conditions for existence cf the 
fundamental frequency forced oscillation. 


A free oscillation frequency response, i.e. F = 0 can be sketched 


from equation (3-23), any value of force function other than zero also can 
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be sketched on both sides of the response curves of F = 0. It is to be 
noted, that the phase angle between the force function and the error signal 
is opposite, i.e., 8 = x, when the response curves are to the left of curve 
for F = 0, and the phase angle between them is in phase, when the response 
curves are to the right of it. In other words, it is according to whether 
the frequency is less or greater than the frequency of free oscillation for 
that particular amplitude of error signal with a constant of force function 
input. 

A typical error vs frequency (EB, 7% cv ) response curves and phase re- 
lation sketch for a different value of forcing function and different char- 
acteristic of nonlinearity are indicated schematically in Fig. 3-5. The 
response curve for free oscillation, corresponding F=0, is drawn as a 
dash line. 

In this respect the behavior of the nonlinear oscillation is the same as 
that of the linear oscillation. One sees that the response curve in the 
nonlinear cases could be thought of as arising from those for the linear case 
by bending the latter to the right for a hard spring saturation nonlinear-~- 
ity, and to the left for a soft spring saturation nonlinearity. 

A response curve in "Input vs Frequency" plane ( e. vw ) can be 
also obtained from relation of: 

6 (t) = 6 (t) - e(t) (3-27) 
or: 


6 (t) # F Cos( wt + 8) - E, Coswt (3-28) 


In the case of @ = 0, equation (3-28) becomes: 
Q(t) = (F - E,) Coswt (3-29) 
Hence: 


9 =F-E (3-20) 
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OT: 


om = F - Q. (3-31) 
Inserting equation (3-31) into equation (3-21): 
2 3 
Ka, (@ -F) += Ka, (@ - F)” 
ame (it (3-32) 
Cc 
For the case of 6 = x, equation (3-28) becomes: 
6 (t) = (F + E,) Cost (3-33) 
Hence: 
- -34 
Ey . (F + 6.) (3-34) 


Inserting equation (3-34) into equation (3-23): 


3 3 
K a, (6. + F) + Z 84 (6. + F) 
we = 2 Ey ES ees Ge a LE We ee a 
6 
Cc 


1/2 (3-35) 


From equations (3-32 and (3-35), if F = 0, the result is the same form 
as equation (3-23), that means the response curve for F = 0 is the same as 


plotted in the E,. UW plane, only the difference is that the value of E 


1 


changes to oe The response curves for the values other than F = 0 are 


] 


changed, and the phase relationship between the input force function and the 
output is reversed, 

There is a numerical example for plotting the response curve in Ew 
and 2. v W plane respectively in section 3-7 of this chapter, and also 
experimental response curves are plotted. 

From the analysis of the above and the numerical example in section 3-7, 
it is shown that a saturation nonlinearity feedback control system can have 
fundamental frequency oscillations. This not only depends on the character- 
istic of the nonlinearity, but also depends onthe amplitude of input force 
function, 

3-4 Forced Oscillations: Saturation Nonlinearity with Damping. 
As explained previously, if there is no damping present in the system, 


there is an oscillation error signal e(t) = E,Cos(@t) either in phase or 


1 
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out of phase of 180 degrees with the forcing function of input. In the 
case where damping is present however, the forced oscillation displacement 
and the forcing input can be expected to be out of phase, just as in the 
corresponding of linear system. 

Recall the basic equation (3-6) and the restoring force function equa- 


tion (2-28): 


An 


E+/E + KE(E) = -A Cos(wt + 6) - BSin(: t + 6) (3-6) 


f(E) = a,E + aE (2-28) 


Then equation (3-6) becomes: 


B+ 28 + Ka,E + Ka ,E =» =AGos(wWt + 6? - BSin(WMe -.0)  °(3=36) 


Recall one solution of equation (3-36) is: 


e(t) = EB, Cosvt (3-15) 


Inserting equation (3-15) into equation (3-36): 


2 3 3 i es 
(= (i Ey + Ka jE, + i Ka Ey yCoswvt ol WE Sin vt + a Ka Ey Cos3c.t 
= (-ACos@ - BSin@) Cosa’ t + (ASin®d -B Cos@) Sin. t (3-37) 


Neglect the higher order terms of harmonic and equates the coefficient 


of Sin .)t and Cost respectively, hence: 


2 | 
(= Ww *E, 4 Ka,E, aw Ka,E, = = ACos@ - BSin®@ (3-38) 
wu By zx A Sin® - B Cos®@ (3-39) 
Squaring equations (3-38) and (3-39), and adding: 
oh 3 S02 ~ 2 7 2 
arn, = 224) Be S 
(= tc Eq, + Ka jE, + Z Ka ,E) ) +x WE, =z=A +8 (3=40) 
or: 
= 3 gee 2 2 ~ eee ye a 
(-WE, + KajsE, +7 Ka,E,) +4 WE, 2 UF te OF (3-41) 
also; 
“, 2 2 cS 2 3 2 zy 2 ae 
yi e — are 4 = ot 
WE, - F) - [2E, (Ka, + 7 Kak, ) Cia 
e,2 3 
+ (Ka, +4 Ka ,E oe = 0 (3-42) 


28 





As a check of equations (3~38), (3-39) and (3-41); first consider 
the response curve in the high frequency and the low frequency ranges, 
the phase relation between the error signal or control signal and the in- 
put forcing function is either zero or 180 degrees. In this case, the 
sine terms of equation (3-38) and (3-39) are zero, the response curve is 
the same as equation (3-20). That means for a damped system if the phase 
angle between the input and the output is either zero or 180 degrees, the 
frequency response is exactly the same as the system with no damping. 

On the other hand, if the value of damping is very small, the re- 
sponse curve will be very close to the sketch of those of Fig. 3-5. The 
only difference when very small damping is present is that the response 
curves are rounded off in the vicinity of curve for F = 0. A sketch of 

them is shown in Fig. 3-6. 

For the phase investigation; recall equation (3-39) and rewrite in 
this form: 

Sin@ - .4 Cos@ sik, /F (3-43) 


Or: wl E 
-] me A - 
Sin(@ - Tan -—~ ) = Flaw + 2)” (3-44) 


Hence, the phase angle between the error signal and the forcing function 


of input: we 
-] << 
1/2 + Tan ae (3~45) 


ma | J 
ae rer 
Examine equation (3-45) in the case of «< = 0, @ is either zero or fn, 


when damping is present the value of @ is a function of @ and E if the 


1? 
values of ~ and F are fixed. It should be noted, there are two values of 
phase shift for one value of error signal; one is for the low frequency 
(w<Wn) and the other is for the high frequency (w >wW,,); one is from 
equation (3-45) and the other ig from the 180 degrees minus the value of 


equation (3-45). 
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For the existence conditions for forced oscillation when damping is 


present, recall the equation (3-42), and put: 


2 
Ka, + + Ka ,E, = X (3-46) 


and: 
ya 


= a -[2B,°x ae ce,” = F*)] + EY X 2 Q (3-47) 


2 
1 


Solve equation (3-47): 


(E 


— 
- 


W == ae i | aT ad (EPA 1: of (eee. iL BAO FE 2G ~FO) (3-48) 
It is to be noted that the conditions for forced oscillations in an 
undamped system are still valid for a damped system within the low and 
high frequencies, in other words, it is valid for the condition of phase 
angle either zero or X. 


From equation (3-48) where forced oscillation exists, the value of 


should be real, hence; the condition for equation (3-48) to be real: 


E, - BS 0 (3-49) 
and; 
2B X = (i - B*) aa) (3-50 
mr 3 
X | 
X > a | (E,” . F*) > 0 (3-51) 
2ZE 
1 
and: 
, 2 2 
[2E 5 Ee 6 3 - F*) | oes GE 247 Op ar F*) >o (3-52) 
1 l | l 1 
or: 
ey. 
4B OX pe 4 st on Aeon 4B, x A (E, - F’) (3-53) 


Equation of (3-49); (3-51) and (3-53) are the conditions for the 


forced oscillation of a damped system, 





3-5 Jump Resonance: 

Recall a sketch of frequency response curve of a soft saturation non- 
linearity shown in Fig. 3-7. ‘Vhen a experimental frequency response is 
being measured for a system, this measurement can be done by an analog 
computer to measure the amplitude of the input and the output, at the same 


time measure the phase relation between them. 
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Fig. 3-7 The Region of Jump Resonance 

It is found that under certain conditions of amplitude and frequency 
an infinitesmal change of either frequency or amplitude of the input signai 
causes a large and discontinuous jump in amplitude of the output, at the 
same time, a discontinuous phase jump occurs. Fig. 3-8 shows a set of 
amplitude jump and phase jump of a typical system from an Electronic Analog 
Computer experimental work. 

Fig. 3-7 is shown as a typical system with a soft saturation non- 
linearity and with a fixed amplitude of the input, which is larger than the 
value of saturation of nonlinearity. If we are going to measure the frequency 
response, being with a low frequency, then increase the frequency, given a 
response from A to B. At point B with a frequency ww, , an infinitesimal 
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increase in frequency causes the amplitude to jump upwards from point 8 to C, 
and with a phase shift from below 90 degrees to a value above 90 degrees. 
Further increases in the frequency cause the output to follow the gradual 
smooth curve to D. Now, if the input frequency is reduced from higher 
values, the output will follow the smooth curve DC. When the frequency is 
reduced to the point E, if the frequency is infinitesimal reduced, the out- 
put jumps down from point E to point F. 

It is found that either from the theory analysis or from the experi- 
mental work, a hard spring saturation nonlinearity can also produce the 
jump resonance, only the difference between them is the jump in the oppo- 
site direction from each other. A typical system for a hard spring satura- 
tion nonlinearity will be shown in the experimental section. 

From the above discussion, the portion of the curve from point B to 
point E cannot be obtained experimentally, usually this portion of the 
curve corresponds to an unstable operating condition. 

3-6 Lower and Upper Jump Frequency of the Jump Resonance: 

As previously discussed, the jump frequency of ‘%, is called "Upper 
Jump Frequency" and the jump frequency of ce’; is called ''Lower Jump Fre- 
quency". Recall the equation (3-41) and the sketch of the response curve 


shown in Fig. 3-7: 


2 
1 


The equation of the locus of the vertical tangents can be found by 


(-W"E, + Ka,E, + 3. KaE, >)” eel” EB =i Ee +l F (3-41) 


differentiating equation (3-41) implicitly with respect to E, and setting 


1 
dw /dE, equal to zero, the result is: 


+ oe 3 2 oe 9 2 2. a 
ae t pea aarti jaa me — 2 / 
(~ uw + Ka, el Ka,E, )¢ + Ka, + 7 Kae, rat = (3-54) 


Check the equation of (3-54), if the damping is very small, the term 
xu ~ can be neglected, that leads to a pair of equations: 


3 3 
= aE, ) = 0 (3-55) 


wo K(a, +7 
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W = Ka, + F a,8,°) = 0 (3-56) 
The curves corresponding to these equations are shown in Fig. 3-9. 
It is to be noted, the equation of (3-55) is the same as the response 
curve equation (3-23) for the free oscillation of an undamped system. The 
equation of (3-56) is the locus of the vertical tangents of curve of Fig. 
B=. 
In the case when damping is present, the locus of vertical tangents 
equation (3-54) will be directed by the equations of (3-55) and (3-56), 
in particular there should be one branch near the response curve for free 
oscillation, and another near the points where the curves for the undamped 
system, In other words, the curve must appear as they are shown in Fig. 
e=10. 


Rearrange the equation (3-54), and solve for « , which leads: 





Oe ee enema ore 


i 
ts ee, = a & r, - _ ; = ate 
Oa | Zend, 3 KK, xe?) x yf CONG: 7A Eh 


Or: 


_ a A. ee 
w= |b Cou Beas tett) 2 f Ceals 2 HtasbeoL)* 40" 30 1GG EEO 3058) 





ee ee are ee me a o 
= — ~ -—me 


In which ... =; Ka,, it is the natural frequency of the linear system. If 


Ee 
we know the constants of system, and the amplitude of fundamental fre- 
quency forced oscillation, the upper and the lower jump frequency can be 


calculated from either one of the equations (3-54), (3-57), and (3-58). 


3-7 Analog Computer Analysis for a Second Order System with Soft Saturation 
Nonlinearity: 


Consider a system with a block diagram shown in Fig. 3-ll. 
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Fig. 3-11 Block Diagram of a 2nd Order System 
with a Soft Saturation Nonlinearity 


In which ne restoring force function of nonlinearity is: 
£(E) = 2E - 0,032E> (3-59) 
The characteristic curve of restoring force function is drawn in Fig. 3-12 
where: 
E# 1,00 1.50 2.00 2.50 3.00 3.50 4.00 5.00 6.00 
feb @wi.97 2.89 3.75 64.50 95515 5.65 5.95 6.00 5.10 
Before the experimental analysis consider the theoretical calculations 
and curves. First consider the response curve on the &. 4 plane for the 
system without damping, and with an input function: 
6 (tc) =z 3 Cos(wt + 6) (3-60) 


Recall the equations of (3-21); (3-22) and (3-23) and leads: 
Daly 2 


w= (4 - 0.0485, ) (3-61) 
for free oscillation and; 
3 
ae 7 2 ae, (3-62) 
for the phase angle 6 = n, and: 
4E, - 0.048E,° 1/2 
cl = ( (3-63) 
SS, | 


for the phase angle 6 = 0. 
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The calculations and curves of equation (3-61) to (3-63) are shown 
in Table 3-1 and Fig. 3-13 respectively. 


Table 3-1: Calculation for response curves on the Ew plane: 


Amplitude of Frequency for Frequency for Frequency for 
Error Signal Free Oscill. @=nx 6 = 0 
(é,) (« ) (a! ) ( #') 
0.00 2.00 / / 
0.50 1,996 0.75 / 
1.00 ARE, 0.99 / 
2.00 La2> 1o23 / 
3.00 1.89 | / ‘ 
3. 5U L.&4 1,35 4.88 
4,00 1.80 1,36 3.60 
5.00 1.68 lee 2.65 
6.00 eye Lec 2586 
7,00 1,29 1.08 1.73 
8.00 1,00 0.80 Loe 
9.00 0.35 0,28 0.42 


For the calculation of response curves and phase shift, when damping is 


present, recall equation (3-42), and leads: 


2E- 2 
: 2 1 2 EL cee 
Kit 13 0.25°— (4 - 0.048E,°)] + zz (4 -0.048E,)~ = 0 
ee li 2 2 l 
Pry -3 EY - 3 
(3-64) 
for the phase shift, recall equation (3-45): 
eal By ud -] 
@ = Sin Flore we il? + Tan iy (3-45) 


It is to be noted, from equation (3-45), and the response curve for an un- 
damped system shown on the =, ~~ plane, when the frequency is larger than 
the frequency of free oscillation with a particular amplitude of input, 
the phase angle ® is zero and is n, when the frequency is below the free 
oscillation frequency. Hence; the phase calculation from equation (3-45) 


should be modified as: 


aie 1 0225 
6 = 180 - (Sin — ee: Tan i (3=65) 
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for jlieculy , and: 


0.,25E ? -|}---— 
-] 1 + San (3266) 


3( 47+ 0.257) 


for wm > és. 

If the response curve is presented on the “ 1. plane, the condition 
of equation (3-65) and (3-66) should be reversed. 

The calculation for the system with damping of 0.25 and the curves are 
shown in Table 3-2 and Fig. 3-14 respectively. 

The jump frequencies of system can be calculated from equation (3-57) 
and also can be seen from the response curve. It is shown schematically 


in the diagram of Fig. 3-14. 


Table 3-2: Calculation for Response Curves (with a damping of 0,25) 


Amplitude of Frequency Phase for Frequency Phase for 

Error Signal for w <wi ta)“ Wig fore >t LD) ay 
(&,) (x) Ce (uw) Cv) 
1,00 1,00 164.00 / / 
2.00 1,24 161.00 i i 
3.00 see 158.50 / / 
3.50 1,36 156.50 4.83 3.00 
4,00 1.37 155.00 3.58 7.00 
4.50 1,40 153.00 2.80 12,00 
5.00 1. 36 151.50 yar 15.50 
6.00 1,24 144,50 2.05 21.00 
7,00 1,09 135,00 1.64 28.50 
8.00 0.87 121.50 1.16 45.50 
8.70 0.65 90.00 0.65 90.00 


For the experimental analysis, the first investigating is the restor- 
ing forcing function, It is to be approximated as a two slope lines, one 
slope is 2, and the other is zero. The circuit setting in the analog com- 
puter is shown below, and also ali components to be used are schematically 
in the circuit of Fig. 3-15. 

The characteristic curve for the restoring forcing function from the 
computer is shown in Fig. 3-16; (a) is shown the waveform of input (e, ) 


and output (e)>» (b) is shown the characteristic of the output to the input. 
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Fig. 3-15 Operational Diagram of a Function 


Generator for f(E) = 2E - 0.032E° 


The operational diagram of analog computer setting up for system of 


block diagram in Fig. 3-11 is shown in Fig. 3-17, in which: 


e(s) = %jM(s) - 5. (8) (3-67) 

ie) - 22 leery (3-68) 
Aw 4,5 rf 

6s) - aa es. (s) (3-69) 


From equations (3-67) to (3-69) system scaling is shown below: 


where: ae = 1,00; x22 1.00; 2281.00; 2.00 
wpa, Re /Ry = 1.00 a, * 1.00; Re, = 1.00; R, = 1.00 
W, = a, Re /R, = 1.00 a, = 1.00; R, = 1.00 
W. = R.AC 
3 £2 £1/4, = 4.00 a, = 0.50; Reo x 1,00 
w a,Reo/a,R, = 4,00 a, = 1.00; R, = 0.50; 
We = ag/ RoC eo =: 2.00 es = 1.00; Re = 0.50 Ceo = 1.00 
We 2 a Re /Re = 1.00 ap = 1.00; Re x 1.00 Re3 sz 1,00 


in which R in megohm; C in //f 
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The results from the analog computer experimental work is condensed 
in Table 3-3; and the response curve with phase shift represented on :, 
plane is shown in Fig. 3-18 respectively. Also the jump frequencies are 
indicated schematically in these diagrams, 

A theoretically calculated response curve in ~ . plane for comparing 
is shown with the experimental response curve of Fig. 3-18. Comparing the 
results from the computer experimental work with the results from theoretical 
calculation, it is shown that the results are very close within high fre- 
quency and low frequency range. There is 4 resonance phenomena between the 


jump frequencies. 
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Table 3-3 Data for Response Curves from Analog Computer for System 
= 0.25; K = 2; £(E) = 2E - 0.032E3: G&(t) = 3,0Cos(.t + @) 


(a) Frequency I[ncreasing: 


Frequency of Amplitude of Phase of Amplitude of Phase of 


Tnput Error Signal Error Out put Output 
(W ) C= 7 (fo ) (@ ) (co) 
0.500 0, 38 180.00 320 0.00 
0,628 0.50 180.00 3.40 0.00 
OER fee ee, 0.75 175.00 3.60 5.00 
0.880 Le25 168.00 4.00 9.00 
1.000 1.50 163.00 4,50 12.00 
Ls £30 2.15 158.50 5.00 24,00 
1,260 3,00 150,00 5.40 35,00 
1,260{+) 9.00 / 6.50 / 
1.380 Sao 25.00 3620 140.00 
Ee LO Teel), 14.00 4,50 155.00 
1,640 7.00 5.00 3.80 159.00 
1.760 6.50 0.00 3.40 165.00 
1.890 6.20 0.00 3.00 172.00 
2.010 6.00 0.00 2.60 180.00 
2.140 5.50 0.00 2.20 180.00 
2.260 5.20 0.00 2.00 180.00 
2.390 00 0.90 1.60 180.00 
2.510 4.80 0.06 1.40 180.00 
2.830 4.40 0.00 1.20 180,00 
3.140 4,00 0.00 1,00 180.00 
3.460 3.80 0.00 0.80 180.00 
3.770 a7 0 0,00 0,60 180.60 
4.090 3.60 0.00 0.50 186.60 


(b) Frequency Decreasing: 


Frequency of Amplitude of Fhase of Amplitude of phase of 


Input Error Signal Error Output Output 
Cw ) (E,) ( i) (G) ( ©) 
250 4.75 0.00 1.50 180.00 
2,200 54:20 0.00 2.00 180.00 
1.890 6.20 0.00 3,00 180.00 
1,760 6.80 0.00 5220 173.00 
1.640 7,10 0.00 3.80 168.00 
1.510 7.58 ,00 4.50 162,00 
1, 380 &. 30 12,00 5.50 156.00 
1,260 9.20 23.00 6,80 140.00 
ie 1 30 10.00 34,00 8.00 135.00 
1,000 11.50 65.00 9.50 120,00 
0.880 11,00 i 8.70 / 
0.880(-) 0.75 175.00 4.00 0,00 
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The resonance phenomena in Fig. 3-18 is actually similar to the re- 
sonance in the linear system, [It depends on the value of damping in 
the system. Another experimental result for the same system and same 
value of input, only changed the damping from 0.25 to 0.5, and is shown 
in Table 3-4 and Fig. 3-19 respectively. In this case the resonance pheno- 
mena is much lower than the case of a damping 0.25. 

It is to be noted, if the damping of a system is increased the range 
of jump frequencies is decreased, For this reason, when the damping in- 
creased to limit value, at which no jump phenomena existed, it means that 
the system becomes a linear system. On the other hand, if the system is an 
undamped or a very low damping system, it will be very oscillatory or a 
limit cycle exists. 

The response curve represented on "Output vs Frequency" (0.7 uw ) 
plane can be calculated from equations (3-32) and (3-35). An experimental 
response curve shown on the ¢-» plane from the system with a damping of 
0.25 is shown in Fig. 3-20. 


3-8 Analog Computer Analysis for a Second Order System with Hard Saturation 
Nonlinearity: 


Consider a system with a block diagram shown in Fig. 3-21: in which 














(9, (5) Ss CS) Scie ef oi "1 } | ? a Oe ($) 
| ie SOS rbizh) 
(. e 





Fig. 3-21 Block diagram of a Second Order System 
with a Hard Saturation Nonlinearity 


assume the restoring forcing function of the nonlinearity is: 


£(F) = 0.425E + 0.005E> (3-70) 
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Table 3-4, Data for Response Curve from Analog Computer for System 
ef = 0.50; K = 2; £(E) = 2E-0.032E3; 6,(t) = 3.0Cos(U!t+e) 


(a) Frequency Increasing: 


Frequency of Amplitude of Fhase of Amplitude of Phase of 
Input Error Signal Error Output Output 
(u ) ( &,) (¢) ( 6.) (¢) 
0.628 0.50 173.00 3.40 7.30 
0.755 0.75 165.00 3.50 12,00 
0.880 1,00 155,00 4.00 16,50 
1.000 1.50 140.00 4.75 24.00 
1.130 Led 125.00 6.00 30.00 
1.260 8.30 91,00 8.00 85.00 
1.380 8.10 47,00 6.25 120,00 
1.510 - 7.50 36.00 Sieo 135.00 
1.640 7.00 25.00 4.25 148.00 
1.760 6.50 20.00 3.50 160,00 
1.890 6.00 12.00 2.75 172.00 
2.010 5.50 5.00 A / 
2.140 5.00 0.00 / i 
2.260 4.50 0,00 . lv 5 180.00 
2.390 Wn 0.00 / / 
2.510 4.00 0.00 1.25 180,00 
2.810 3.75 0.00 1,00 180.00 
3.140 3.50 0.00 0.75 180.00 
3.460 So a 0.00 0.50 180.00 
on 77/0 3.25 0,00 0,30 180,00 
4.090 3.10 0.00 / / 


(b) Frequency Decreasing 


Frequency of Amplitude of Phase of Amplitude of Phase af 
Input: Error Signal Error Output Output 
(w) ( <,) (2) (6) Co) 

2.510 4,30 0.00 Ve2o 180.00 
2.200 5.00 0.00 NO. 180.00 
1.890 5.75 0.00 2.75 172.00 
1,760 6,25 7.50 3.50 160.00 
1,640 7.00 15.60 4,10 148.00 
210 7.75 27.00 5.00 132,00 
1. 380 Seo 32.00 6.50 120.00 
1,260 8.50 60.00 7.50 97.00 
1,130 7.00 80.00 8.20 40.00 
1.000 2.00 130.00 4.75 19.50 
0.880 1,25 145.00 4.00 0.00 
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The characteristic curve of equation (3-70) is shown in Fig. 3-2la, where: 

Bae 1,00; 2.00; 3.00; 4.00; 5.00; 6.00; “sauce 8.00. 

F(Ey 00,43; 0.89; 1642; 2502; 2.7534 “Se03; 4568") S26. 

A theoretical analysis can be calculated from equations (3-70) and 
(3-21) to (3-23), as we did in section 3-7, Only the analog computer anal- 
ysis is considered in this section. 

For the simulation of the restoring forcing function of nonlinearity 
from the analog computer, it is approximated by two straight lines, on 
with a slope of 0.425, and the other with a slope of 2,00. The circuit set- 


ting in the computer is shown in Fig. 3-22, and all components to be used 


are schematically in the diagram. 


y —, 


: ka A tt |. | = ~ Ae le | 







Ps tas 


i, | 


Fig. 3-22 Operational diagram of the Function Generator 
for f£(E) = 0.425E + 0.005E2 


The characteristic curve of restoring forcing function from the com- 
puter is shown in Fig. 3-23. 

The operational diagram of analog computer set up for the whole 
system is mostly the same with the diagram of Fig. 3-17, the only dif- 


ference is the function generator circuit. 
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The results from the analog computer are condensed in Table 3-5, 
and the response curves with phase shift shown on the ¢;-* and - » 7 
plane are shown in Fig. 3-24 and 3-25 respectively. In which the input 


forcing function is 8 (t) = 3.50 Cos(c t + @). 
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Table 3-5. Data for Resporse Curves from Computer for a System with a 
Hard Saturation Nonlinearity, in which: ~ = 0.25; K = 2; 
f£(E) = 0.425E + 0.005E3; Q(t) = 3y5°Ces (irte4-8). 


(a) Frequency Increasing: 


Frequency of Amplitude of Phase of Amplitude of Phase of 


Input Error Signal Error Output Out put 
(Ww) Cay) (O ) ( ©) (@) 
0.500 1.75 180.00 5.00 0.00 
0.628 35 30 180.00 6.50 4.00 
0.754 6.50 175.00 9.00 7.50 
0.880 8.00 168.00 10.50 10.00 
1.000 9.00 164.00 12.00 11.50 
1.130 10.50 158.50 13.00 13.00 
1.260" 2.25 152.0 14.50 15,00 
1.380 14.00 145.50 16, 30 17.00 
1.510 16.25 141.00 18.75 20.00 
1,640 19.00 134.00 21,30 24.00 
1.760 eee 124.00 22,00 37.00 
1.890 26,00 102.00 25.00 45.00 
1.950 26.50 / 25550 55.00 
2.010 26.00 24.00 25.00 90.00 
2.140 4.50 20.00 0.86 150.00 
2.260 4.35 15.00 0.65 165.00 
Z.520 4,25 7.50 0.50 180,00 
2.830 4415 0,00 0.40 180.00 
3.140 4.00 0.00 4.30 180 SO - 
2.770 3.80 0.00 / i! 


(b) Frequency Decreasing: 


Frequency of Amplitude of Phase of Amplitude of Phase of 


Input Error Signal Error Output Output 
(i ) ( 5) (6 ) ( ©) ( «€) 
2.310 face 5.00 0.50 180.00 
2.400 4,35 ae 0.63 180.00 
2.270 4.50 13.50 0.68 172.59 
2.140 4,60 15.00 0.75 170,00 
1.890 4.85 21.00 | os 157.50 
1.760 5.00 23.50 125 155.00 
1.640 5.35 25.90 1.60 152.00 
210 AT 26.00 eS 150.00 
1.380 6.75 27.50 16,30 18.50 
1.320 ees / / / 
1,260 110 150.00 14.50 15.50 
30 10.25 155.00 13.00 13.50 
1,000 9,90 160.60 11.70 12500 
0.880 7.50 172,00 10.50 10.00 
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CHAFTER IV 

SUB-HARMONIC OSCILLATION IN NONLINEAR FEEDBACK CONTROL SYSTEM 
4-1 General Description: 

Up to the preceding chapter, there is only discussed the forced 
oscillations for which the frequency is the same as that of the external 
force. A nonlinear control system, under suitable conditions may have a 
steady state oscillation in which the main component has a frequency which 
is a fraction 1/n (n is a real and positive integer) of the forcing func- 
tion frequency. These oscillations are called '"Sub-Harmonic Oscillation". 

Recall a linear system (Section 2-3), there are two oscillation terms 
in the transient response, if the frequency of the free oscillation is 
( /n (n is an integer), then a force function of frequency ~ can excite 
the free oscillation in addition tothe forced oscillation of frequency « . 
But as t increases, e ti t/2 damps the free oscillation until only the 
steady state oscillation remains. Hence there is no sub-harmonic or super- 
harmonic oscillation existing in a linear control system at steady state. 

In the case of nonlinear control system, the steady state response 
may not have the same frequency of oscillation as the frequency of the 
forcing function. In fact, the frequency of the system response depends 
not only upon the frequency of forcing function but also unon the amplitude 
of the input. The principle of superposition no longer holds and the linear 
circuit characteristic is not applicable. In general, it can be said that 
the nonlinear system can have wide variety of almost periodic oscillations, 
the frequencies of which differ from the frequency of the applied force and 
vary with time as well as with different initial conditions. 

Sub-harmonic oscillation always existed with a frequency higher than 


the frequency of free oscillation with a particular amplitude of error signal. 
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The response curve which is a part of an ellipse or a hyperbola in the 
Ey plane is always nearly parallel to the right side of the free oscil- 


lation response curve. A typical sketch of a soft saturation nonlinearity 


is shown in Fig. 4-1. 
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Fig. 4-1 Fundamental and Sub-harmonic Response 
Curves for a Soft Saturation Nonlinearity 


Under some suitable conditions, there may exist any order of sub- 
harmonic oscillation, or only one type of order can exist depending on 
the characteristic of the nonlinearity and the input forcing function. 

Sub harmonic oscillation and response of nonlinear vibratory systems 
are often discussed in the field of mechanical system as There is no 
body of knowledge applied to feedback control servomechanism. For the feed= 
back control system, there is a standard differential equation (3-6a), it 
is very similar to the equation for the mechanical system, the only dif- 
ference is that the forcing function of a feedback control system is a func- 
tion of input frequency. 


Although any order of sub-harmonic oscillation can exist in a non- 


linear feedback control system, in this paper only the sub-harmonic of order 
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2 and 3 are discussed. 
4-2 Sub-harmonic oscillation of Order 2 without Damping: 
Sub-harmonic oscillation of order 2 can exist in most nonlinear feed- 
back controls, When it occurs, the frequency of output (¢€.) or error 
(2) ts only 1/2 the frequency of input ( “-). A typical diagram from 
the computer experimental work shown on the "Input vs Output" plane and the 
waveforms of "Input and Output" are shown in Fig. 4-2 and 4-3, respectively. 
Recall the basic equation for forced oscillation in nonlinear system 
with un-damping Eq. (3-11). 
E + KE(E) = -A Cos(& t + 0) (3-11) 
For the same type of restoring force function with chapter 3, leads: 
Et Kae + Ka,E” = - A Cos(@t + 6) (3-12) 
Assume one solution of equation (3-12) is: 
Roe Cos swe + Ey Coswt (4-1) 


Where E is the amplitude of sub-harmonic of order 2, EB is the amplitude 


1/2 


of fundamental frequency, since: 


— 


i «ee > 
E = Zé E4/2 Cos zwt - wv Ey Coswt (4-2) 
a ee a ey eee eal! 
- = & ly, a ae - 4 iy :, Jie ae > (NS -. i + 2F JCOSWE 
a ee ; 3 . 
ta (by a A ee eS a a Ky 
_ — bee 
a te } : = me: =, Cia 3:98 (4-3) 


We are only interested to the terms of second sub-harmonic and fund- 
emental freauency, inserting equations (4-1) to (4-3) into equation (3-12) 


the result is: 


* " ps 4 4 < - J j aes ied 
a Ey, + Ky TE Kisby, t Eads Ey SY Was We + (- WE, the, 


ft. 


Sed 


f 
iN 


3 a eee 2 cm" ~~ ail Aet 
=> Hy {\ (Oty i, T ’ j y/ sip” = ; i re a. fs ” Z las ; 5 (4-4) 
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As discussed in Chapter 3, Section 3-3, if the frequency is larger 
then the frequency of free oscillation, the phase angle @ is zero on the 


“Error vs Frequency” plane, hence equation (4-4) leads: 


ama! ; w 1 Le pels 
(UY = 4A8, ¢ Bea + b KE, (4-3) 


Ww = i ie RS 
= (4-9) 


Equation (4-5) is the condition for existence of the sub-harmonic of order 


2. Solve Ey/2' 


= — 4 ca Be = al Lf . 4 ie Vs ~ 
Ey = 4 ( 3Kas 3 0; * LES) Ca) 


Since E/2 should be real; which leads: 


WE (Aa, + 18 1s) és v - 
For the first iteration method, consider a. equal to zero: 
Equations (4-5) and (4-6) leads: 
Wr = Ake, = 4h), (4-9) 
B= $1 (4-10) 


Inserting equations (4-9) and (4-10) into equation (4-8): 


WZ 2 ty + B ifay F7)*F dy ZO (4-11) 
From Equation ("-5), it represents an ellipse or a phperbola in the Eys2 
vs « plane depending on the sign of a Ans also «) has a maximum when a. 
is less than zero, and a minimum when a4 is greater than zero. Actually, 
the sub-harmonic oscillation only existed within an interval of frequency. 
That means only one part of ellipse or hyperbola can exist, as shown in Fig. 
(4-4). 
4-3 Sub-harmonic Oscillation of Order 2 with Damping; 

Recall the basic equation (3-41) for a damped systen: 
3 


E+¢dE + Ka,E + Ka,E = - ACos(wt +6) - BSin(wet + 8) (3-41) 
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Assume one solution of eaustiou (3-41) is: 


= =i as Se 7 yA Caswt * ae 7 Lee (4-12) 
Since: 
ey a ha a) (4-13) 
@: ee ‘ adh Tits fy Vedic 1s L7(5 aie t/l 
a Zt Ey, la Zale «Sy Wr Cawe — Byawuwet (4-14) 
ge oy Va a. i te = i = a bs (4 ca t 
= © (gE y tT ene 8 ere 7 ore eh 
eo Yo _ 
f (pe z Lia a ey es pac Gwe 
Bis eZeyn rey Gy) re Wt 
Bae BT me as 15) ¢ (4-15) 


Inserting equations (4-12) to (4-15) into equation (3-41) and equal- 


e l co 
ing the coefficiants of Cos G4E Cosu¢t, and Sinac, give a result of 
we tae ae. UE 
- fH Lyf, Tha, fy % ‘ 5 (fey + Ztm r see ) == U (4-16) 
BE 4 ~ : os — a 
“EaW t ht W ORR Ey oh Os (ZEA t Zip ae a By, 
ooo Ae iy! bs Jas iy (4-17) 
= a sb _—e ye “ 2 Ags 
=Br, — ea A + Kite a, " ARG (A “AT BOE Seg ips? 
=== Ahm! ~ Bce+ © (4-18) 
l 
= (apc, 


From equation (4-18), it is the coefficient of the term of Cos 2 


where Ey/2 #0, hence: 
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=a - , , (4-19) 
Solve for E12? 
Bisa = y ‘ie ? Aig 4 4 ! 4 tT ° ) (4-20) 


The value of Ea y2 should be real therefore, the existence of 2nd order 


sub-harmonic with damping: 


a -_——— 


f 


we ithe, + biw, Sa + baa te a, & 0 (4-21) 





From equations (4-16) to (4-18) for the first iteration method, for 


&. equa | to zero: hence: 

eZ 4ka, (4-22) 
(Ka, - a” JEW, + Kt Ex =z - . FCos® - <« FSin@ (4-23) 
(Ka, - U) Ey, glee Ean = @ FSin@ - «<u FCos@ (4-24) 

Solve ETA and Eis from equations of (4-22) to (4-24): 

= £7 _~_ + — ot een LY s 
ELA S rite Wi rn —)Goh + yee i / (4=25) 
Bo ee LL ef Be ete Dae | (4-26) 
1B one rT = od — 


Inserting equations (4-25) and (4-26) into equation (4-21). A 
general equation for existence of 2nd order sub-harmonic oscillation 
with damping: 


> * 
er ¢ 
an FA SE “4 


*s 
Z 
~-% pe 





fe 2 t; hee r - 4 ry 7 
y= ae cs ge , — ; rg. ‘a i ‘ , 
baa 2 | a7 a 2 Kis fF rh (/ pos ee _— Gs 704 + fg 4 4) 
’ % 2 , ae 
: ’ a cy ey 7 ew 
+ oKagp yy — aoe ee Cee sem eam eye | a 
We a A 4 (4-27) 


For the case of @ equal to 0 or n, i.e., Sin@ equal to zero equation 


(4-27) becomes: 


n= ny A —- a? Ve 
UW ee we ’ i i J — A (3 } a7 si aia f . | i C kas “a 4 1g- (4-28) 
: # ae tt . 
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From equations (4-11) and (4-28) they were shown that either the 
system with damping or without it the second order sub-harmonic oscil la~ 
tion only existed with a frequency less than the frequency of two times 


the natural frequency of linear system for the case of a, less than zero, 


3 


and greater than this frequency for the case of a, greater than zero. 


5 
4-4 Sub-hamnonic Oscillations of Order 3 without Damping: 

Typical waveforms with jump phenomene of 3rd order sub-harmonic 
forced oscillation is shown in Fig. 4-5; (a) shows the transient from 
the 2nd order to the 3rd order sub-harmonic; (b) is shown reversing. 


For investigating the existence conditions of 3rd sub-barmonic forc- 


ed oscillation recall the basic equation of (3-11); and leads: 


E + Ka FE + Ka E> «- A Cos(tt + 6) (3-12) 


Assume one solution of equation (3-12) is: 


] 
= = ie ous 
E E4/3 Cos 3 t + EB Cos&t (4-29) 
in which E13 is the amplitude of 3rd order sub-harmonic oscillation, 
hence: 
E = ee Cos ae t -w E, Cosat (4-30) 
9 “1/3 3 [eral 
3 3 
E = = 2 2 eye 
1 3 2 ? 
re . ev e@@ 06 ad ) 
i (2173 + 6E 1/3 Ey + 3E, )} Cosdit + (4-31) 


From equation (4-29) to (4-31) and equation (3-12) leads: 


( KG ; 3 ) fot 3 - - Nig: Ly { ~ ; a J ra (4-32) 
acy 4 + & a A re r & & — } : } (4-33) 
Equation (4-32) E173 is not equal to zero, hence: 
a + Hebe Tiltenie ~ * = ei (4-34) 
Solve E 1/3 from equation (3-34). 
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o, ad a mm! (4-35) 


The value of E must be real, then: 


3 


Aq : ' 27 hdA3 = ’ 
a Gy 0 ee EL, 7, = ¢ (4-36) 


Equation (4-36) is the condition for the existenceof sub-harmonic 
oscillation of order ? of a without damping system: 
From equations (4-32) and (4-33), by the first iteration method, 


for a. «= 0; then: 


to" = JY Ka, (4-37) 


(Ka, - w )E, = -sF (4-38) 
Solve Ey from equation (4-37) and (4-38): 
9 


Inserting equation (4-39) in equation (4-36) the condition for a 3rd 


sub-harmonic oscillation becomes: 


ae 


Also the response equation of equation 4-3 becomes: 


~N 


oa a ‘| ae =) on 


pe oo eeye 1-ee,.* SS ; =¢ (4-41) 
The equation of (4-41) represents an ellipse or a hyperbola in 


the E173 vs 4 plane depending on the sign of a5, also has a maximum 


when ay 0, and a minimum when a, > 0. A sketch of a 3rd sub- 


harmonic forced oscillation ts shown in Fig. 4-4, 
4-5 Sub-harmonic Oscillation of Order 3 with Damping: 


Recall equation (3-41): 


BE + <b + KajE + Kak == ACos(wt +6) - BSin(wet + @) (3-41) 


Assume one solution of equation (3-41) is: 


ee . : 
E = E1/3008 5 ut + E,,coswt + E ypoiny t (4-42) 
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Since: 


. ] y 


ESE) = Sin guste - BE), & Sinest +E), « Cos mt (4-45) 
E = - Se Cos 5 tL - E,,@ Coseet - - Eo Sinmt (4-44) 
gE = 2 OE + EL/3P 1A 4 2B + 2B, )Cos 5 wt - 

Ze 175 EypSing 4 — V/A (8, 75° Ey 75 Ey + Bees 

3B Ep) Costit +2 Dey ORG + 2 Bi. + Bia) Sinwt + ..(4-45) 


4 
From equation (4-42) to (4-44) and (3-41), by equaling the coefficients 


eh. ; 
of Cosy UC, Sin» t, Cos. t, leads: 


? fq, ~ 4 : 
CM, 12) +2 eee Spee Se 1 ABD (4-46) 
Che, ll .~ r 74 he om és : ; on = 7, + Sr fim ) rele 
. (4-47) 
= 7 ad | "ea. f -— a 
ih & ). . = | : a 4 / 
(AX, vv it “ mM mainte «6OTlUCOU Uf A f , <~" / ‘ rt Ay i - fs ; (4-48) 
y / 7 ' 
Rearrange equation (4-46) and solve it for E173 in terms of ; 
E and E,.? 
1A’ 1B 
é i] ate oa 
~~ / ‘ 9 r We - rtf 
Lf, “eee 2)! etyA TE 7, Be “sey Ke oe (4 49) 
here E/3 should be real, hence; 
P = _ ware @ing rl 
io eFC he, ot eae ES ey (4-56) 
or: 
7 > —— 3 cd /) ite = .—_=- 
ra = 3 y bur i iat ‘a , T ws (4-51) 


Equation (4-51) is the condition of 3rd sub-harmonic oscillation 
existing in saturation nonlinearity with damping system, 


By the iteration method, equation (4-46) to (4-48) leads: 


uf = Ika, (4-53) 
- 4 es - |. i = i. = 

(Ka, w Eva + cu. Eis « FCos 9 = 4! FSin® (4-54) 

(Ka) = (A JELp ~ ji EA = (| FSin @ - /. FCos@ (4~55) 
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Solve equation (4-53) to (4-55), 


EIA a 4 , / os i + = it 7" (4-56) 
Ss es ; . N (4-57) 


_e tf 


Substituting equations (4-56) (4-57) into equation (4-52) and be- 
comes in a general form: 


cel ———— = 





=~ Be JS / aa , a 4 
Pee? pm 7 Ste / 7 ae i | 
uy - | ' oy, saa = a Ley r ce ré. 
——— a 2 ' . ; a) 
ral Oe nO ue WG 7 46 “Bes 
aoe [ot rl — en fa | (4-58) 


In the case of @ equal to 0 or a, i.e. Sin@ equal to zero equation (4-58) 


becomes: 


as = Fi 
4 4 a 
2 ot? biel Zt x 


zx we oa  ' [5.58 F ( /@ pat me a pias 19 i 


~~ 


(4-59) 

In conclusien, the condition for existence of sub-harmonic oscilla- 
tion in a nonlinear feedback control system depends on the characterics 
of the nonlinear element and the amplitude of input. The response curve 
of any order sub-harmonic oscillation with a saturation nonlinearity is 
either an ellipse or a hyperbola which depends on the characteristic of 
nonlinearity. Also the frequency of a sub-~harmonic oscillation existed 
may be either a maximum when a soft saturation nonlinearity is presented 
or a minimum when a hard saturation nonlinearity is presented, The limit 
value of minimum or maximum is almost equal to the number of sub-harmonic 
order times the natural frequency of linear system. 


4-6 Analog Computer Analysis for a 2nd Order Sub-harmonic Osciilation 
System: 


Consider a system with a block diagram shown in Fig. 4-6: 
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Fig. 4-6 Block Diagram of a 2nd Sub-harmonic 
Oscillation System, 


Assume the restoring force function of nonlinearity is: 


£(E) = 2E - 0,0048E° (4-59) 


The characteristic curve for equation (4-59) is shown in Fig. 4-7, in 
which: 

Bam .00;> 2,50; 5,00; 7.50: 10,002. W750; 1aee0. 

fcnpm= 1.995; 4.295; 9.400; 12.97; 15.20; 15.40; 145603 

For the simulation of restoring force function of nonlinearity from 
computer, it is approximated from two straight lines, one of which with 
a slope of 2; and the other with a slope of zero. The circuit setup in the 
computer is identical with the circuit shown in Fig. 3-15. The only dif- 
ference is the value of bias voltage to be used. The characteristic curve 
of restoring force function from the computer is shown in Fig. 4-8. 

An operational diagram setup in the analog computer is almost the 
same with the diagram of Fig. 3-20, only one difference is that the value 
of constant to be used. 

The results from the computer are condensed in Table 4-1 and Fig. 4-9. 
Only the response curve is of interest, hence, only one response curve in 


the error vs frequency plane is presented. It is shown that there are two 
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jumps, one is the fundamental frequency jump, and the other is the second 
order sub-harmonic frequency jump. These jump phenomenas are shown sche- 


matically in the diagram of Fig. 4-9. 


Table 4-1: Data for a 2nd Order Sub-harmonic Oscillation Response 
Curve, in which: f£(E) = 2E - 0.0048E7; K = 4,-/ = 0.25; 
and Q(t) = 2.5 Cos( t + 8) 


(a) Frequency Increasing: 


Frequency of Amplitude of Amplitude of Remarks 

Input lst Harmonic 2nd Sub-harmonic 

(w ) (c,) ( &4) 

0..62& Oxze2 / 

0.754 OF: / 

0.880 0.45 / 

1.000 0.50 / 

1,130 0,63 / 

1,260 1.00 / 

1.380 1, i / 

1.510 LeU / 

1.640 2.00 / 

1.760 2.20 / 

1.890 3,00 / 

1.950 3.25 i 

1.980* 10.50 / *1st harmonic 
2.010 10.00 / upper jump 
2.200 9.00 / 

249.10 7.30 / 

2.810 6,00 / 

2.950 / / 

3.140 3.00 / 

3.460 4.50 / 

3.770 4,00 / 

4.090 cry fe) / 

4.400 3.50 / 

4.710 3a2) / 

5.020 / 4.75% “2ad sub-harmonic 
5.150 / 4.50 upper jump 
5.270 / 3.50 

5.340 2.50 / 
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(b) Frequency Decreasing: 


Frequency of Amplitude of Amp litude of Remarks 





Ce, 


Input lst Harmonic 2nd Sub-harmontc 

(& ) (a, ) ( im) 

2. 340 2.50 / 

5.270 / 3,20 

>.150 / 4.50 

5.020 / 4.75 

4,710 / 5.50 

4,400 if 6.50 

4,090 / 7.89 

3.770 / 9.50 

3.460 / 11.00 

3.140 / 13.50 

2.950 / 14.50% *2nd sub-harmonic 
2.810 6.00 / lower jump 
2.510 P20 / 

2.200 8.75 / 

2.010 10.60 [ 

1.890 1).00 / 

1.760 l2,25 / 

1.640 13,00 / 

1.510 15.00 i 

1, 380 17.00 / 

1,260 21.20 / 

1.130 24.00 / 

1,060 25.50 / 

0.880 0,45* / *1lst Harmonic 
0.754 0.33 f lower jump 
0.628 O22 / 











4-7 Analog Computer Analysis for a ird Order Sub-harmonic Oscillation 
system: 


Consider a system with a block diagram shown in Fig. 4-10: 








Fig. 4-10 Block Diagram of a 3rd Sub-harmonic 
Oscillation System 


in which assume the nonlinearity of system is an "ON-OFF" element, a 
sketch of the characteristic for the computer simulation is shown in 
Fig. 4-ll, it is approximated from three straight Jines; one of which 
with a slope of 15.0; one of which with a slope of 0.1; and the other 
is zero; 

The circuit setting up in the computer is shown in Fig. 4-12, all 
components to be used are schematically indicated with the diagram. The 


characteristic curve from the computer is shown in Fig. 4-13, 
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Fig, 4-12 Function Generator Diagram for an 
"ONe©OFF" Nonlinearity 
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An operational diagram tor setting up the analog computer is shown 
in Fig. 4-14, in which the value of all components to be used are 
schematically indicated in the diagram. The scaling factor for the 
computer setup is the same as diagram of Fig. 3-20 to be used. 

The results from the computer are condensed in Table 4-2, and with 
a response curve shown in Fig. 4-15. It is shown that, there is a pair 
of jumps for each harmonic oscillation response curve, these jump pheno- 
menas within frequency range is accompanied shown with the dtagram of 


Fig. 4-15. 


84 





2 
; 
: 





STI UT UI2eTY YOOTe Fo 
woysfc TOT Gn BuTayZog Toynoneg 
Pie, [RuOlVereds 7 y ote 


3 








= ee 3s NS i —— 









TOU CIOS: 


UOTIOUN 


eee, caer 
a | 
| a “oT Pe! eal | 
| ap) con on) 


j 
} 





85 





Table 4-2; Data for a 3rd Order Sub-harmonic Oscillation Response 
Curve from an "ON-OFF" Nonlinearity System and; .2 = 0.4; 
K = 10.0; 6,(t) = 8.0 Cos(;vt + 8) 


(a) Frequency Increasing: 


Frequency of Amplitude of Amplitdue of Amplitude of Remarks 
Input lst Harmonic 2nd Sub-har. 3rd Sub-har. 
(w ) ( ¢,) ( cy) (zy) 


0.628 0.60 
0.880 leo 
1.130 2225 
1, 380 3.50 
1,640 4.75 
1.760 5.75 
1,890 7.50 
2,010* 50.00 
2.140 46.00 
2.270 42.50 


*1lst Harmonic 
Upper Jump 


2.390 
2.510 
2.630 
2.760 
2.890 
3,010 
3.140 
3.260 
4.270 


39.00 
35.00 
31.50 
28.00 
26.00 
22.50 
20.00 
/* 
/* 


*kTransient and 


unstable 


“Transient trom 
lst Harmonic 
to 2nd Sub-har. 


4,390° 
4.520 
4.650 
4.770 
4,990° 
5.020* 


™~™ ™ ™ ™ Me MM. Me™. MM MM MM MM “MR ME SL 


5.140 
5.270 
5.400 
5.520 
5.650 
ae 70 
5.900 
6.030 
6.150 
6.280 
6.400 
6.530 
6.650 
6.780 


6.910* 


7.030 
7.160 
7.280 
7.410 
7.530 


20.00 
18.00 
16.70 
16.00 
15.00 
13.80 
13.50 
W270 
12,30 
11.90 
11.30 
11.00 
10.60 
10,20 
10.00 


“~~ “™., MM ™ 
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me ™ ™. ™ MM MR MR MR MR MR MM MM MM Mn MR MR MR MR MR Me MR TR TR TM TE, SM, 


35.00 
32.50 
28.50 
27620 
26.00 


*2nd Sub<-har. 


Upper Jump 


*3rd Sub-t. 
Upper Jump 


vr. 





(b) Frequency Decreasing: 


Frequency of Amplitude of Amplitude of Amplitude of Pemarks 
Input Ist Harmonic 2nd Sub-har, 3rd Sub-har. 


Cw ) 


1530 
7.410 
7.280 
7.160 
7.030 
6.910 
6.780 
6.650 
6.530 
6.400 
6.280 
6.150 
6.030 
5.900 
5.770 


5.650% 


3220 
5.400 
5.270 
5.140 
aee20 
4.900 
4.770 


. 4,650 


4,520 
4,390 
4.27C 
4.140 


3.260° 


3.140 
3.010 
2.890 
2.760 
2030 
2.510 
2.290 
2.270 
2.140 
2.010 
1.890 
1.760 


1.640% 


1.510 
1.380 
E50 
0.880 
0.628 


( éi) 


—™ “sy Ne TM, Me, MA Me Te, MA TM. MT, MR MR, TM Me MR. MR, MT, Me MR, ML MRL MR, MRL MM, MM, 


20.00 
22.50 
26.00 
27.50 
31,00 
35.00 
39.00 
42.00 
45.50 
51,00 
35.00 
ay ree 9) 
5.00 
4,25 
3.50 
“30 
1,25 
0.60 


C ®,) 


13.75 
15.00 
16.00 
16.70 
17.50 
19,50 
21.50 
23.00 
25.00 
32.00 
36.00 
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26.00 
27.20 
28.50 
30,50 
32.50 
35,00 
36.50 
37.50 
40.00 
42.50 
43.75 
47.00 
50.00 
52490 
56.00 


m™ TA MA SM, TS TE TS, TE TL TS, MR TL MR, TR TR, TR, STM TM, MR ME, TS MM MM ME SS, 


*3rd Sub-har. 
Lower Jump 


°*2nd Sub<-har. 
Lower Jump 
and Transient 
Unstable Region 


*1ist Harmonic 
Lower Jump 
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CHAPTER V 
FORCED OSCILLATION FOR THE HIGHER ORDER 
FEEDBACK CONTROL SYSTEM 
a-l General Description: 
Up to the Chapter 4, only forced oscillations in the case of the 2nd 
Order system are discussed. [In the extension to higher order systems, 


consider a block diagram shown in Fig. 5-1: 








Fig. S-l Block Diagram of a Feedback 
Control System 


in which G(s) is an equivalent transfer function of the linear element 


of the system, it can be written in the form: 


G(s) = x Qs) 


F(s) (5-1) 
and with a result for the system: 
P(s) E+ Q(s) Kf(E) = P(s) (5-2) 
Where F(s) may be any order of s; as: 
n n-l n=2 ; 
P(s)'= s + bs + bs deieveie ate db, (5=3) 


Similarily Q(s) may be 
m m=] m=2 
Q(s) =s + c¢)8 + CoS + eecccecsoet C. (5-4) 
In general the frequency responsive element 3(s) following the nen- 


linearity is of such a nature as to attenuate the high frequency, that 
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means the order of P{s) should be higher than the order of Q(s). The 
order of system is n. Inserting equations (5-3) and (5-4) into (5-2) and: 


Cg bene r iy “all T ae Wit) 


dee 7 Fi etm 


dea 
t+ & oes, ria ) (5-5) 


Me } 
é é t 


( 
Equation (5-5) is a general differential equation forn of the 
forced oscillation of any order feedback control system, 


5-2 Harmonic Linearization and Iteration Method for Investigating the 
Equivalent Gain of a Nonlinearity: 


“After a restoring force function of the nonlinear element has been 
defined, the equivalent complex gain, N(E,) of the nonlinear element is 
determined by the ratio of the complex amplitude of the first harmonic 
of the output to the amplitude of the harmonically varying input quantity. 

Recall the block diagram of Fig. 5-1, if the system is a 2nd order 


system and with a restoring forcing function; f£(E) = a,E + aE , and: 


K NCE) E 2 @ 
os (8s +H) 7 | So 
or: 
E +.2E + KN(E)E=6, +0, (5-7) 


Assume the input function is 
V, = F Cos(u t + 6) (5-8) 
and one solution of E:; 
E=E, Coswt (5-9) 
and with a result; 


(io-E, +K N(E,)E,)Coswt + / E,Sinwt = -A Cos(ut + 8) 


1 
~ BSin(wt + 9) (5-10) 


Comparing the result with the equation of (3-42) the equivalent gain: 


3 
N(E,) = a, +7 ay (5-11) 


1 





For a general torm expression; the equivalent gain of a non- 
linearity in terms of the amplitude of input ts develrped by Ya. Z, 


Tsypkin. The resuJt for a symmetrical nonlinearity with input of equation 





(5-9) is: 
es 
= . ( - 
N(E,) 3E, Hi £(Z,) + goo (5-12) 
If more accuracy is required, the equivalent gain N(E,) may be: 
ae l i , = 
N(E,) "32, L £(E,) + E(B ,/2) + 3 €C 73 B,/2) (5~13) 


For the case of a nonlinearity with a restoring force function, 
f (E) : a,E + ae, the results from equations (5-12) and (5-13) are ex- 
actly the same as equation (5-11). Hence, for any nonlinearity the case 
of first harmonic input the equivalent gain can be expressed either in the 
form of equation (5-12) or (5-13). 

Equation (5-12) or (5-13) not only substantially simplifies the 
calculation for a second order forced oscillation, it is also probably 
a most important for investigating self oscillation or forced oscillations 
in a higher order system with single or multiple nonlinearities in a 
feedback control system, 
5-3 Forced Oscillations for a Higher Order Feedback Control System: 

Consider a system with a linear transfer function: 

G(s) = OC (5-14) 
s(s + p,)(s + py) (s + Pp.) 


or in the form: 


K(s + c¢ % 
cms ee (5-15) 
Ss + bs + bos + b,8 


From the equation (5-5), the differential equation for the system 


of equation (5-15): 











4 3 2 ea ee oe 
CaBigper mea + byes + by eS ae * I 
dt dt dt: 
4 3 2 
d J d d_) 6 
* ( — 4). =) ee ee , 
male Maat 2 dt 3 dt (5-16) 





in which, 


>) 
{ 


<3 E Cos ue Ct 


) = FCos(ut + 6) 


and; 


Pare 3 
f(E) = a, E ,Cos vt + 7 ane, (3Coswt + Cos3. t) 


Inserting these equations into equation (5-16) and with a result: 


P , | in 
(FUT - bev rid rg ha,c ce * J, ee ale 


, h 43 i? O/, din Qt 


r (0, 8,0°-bhj ec al—-het,£. tv 


7 ; 
f “5 K &@; ee £ Oe 34/6 ~ - Wa, WM Rin 340 


; 5 : - \ os 
= Fut w (ud ~-he dios? +( bt —- 6; )am8 [ORM 


-Palsa (Ww -bs hw —(b6,W0— 6.) CeO) dm 
(5-17) 


Neglect the higher order harmonic terms and equate the coefficients 
of Coswt and Sinit: 


od ‘ > - =~ 
fae - bh CRP RAE 1 = has & 


wef 


it { ( bo bP a ene, Cb. ie bs dink | (5-18) 


» 


E,(06;W*-b,) —(na, + 2- Ka aay 


~bu 


Pel (eb. jweaf —- Ch Wb) eso] (32179) 


Squaring equations (5-18) and (5-19) and adding: 


> : ina 4 
~w ; - 74, = Pr ee i ! / e & WF j 
ECA bla + Kaye to NA 2, EW [CBW 00, Aha; 


= pip] wre re)” -r Cb, tv- k,) J 


a 


(5-20) 
Equation (5-20) is the equation of response curve of the system of 
equation (5-15). If the constants of system are know, a response curve 
in the <,vé) plane can be plotted. 
For the investigation of conditions for existence of forced oscil- 


lation rearrange the equation (5-20) and put: 
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Hence: 


OT; 


The condition for existence of forced oscillation is found frem 
equation (5-22) by means that the solution of «© should be real. It is 
to be noted that the value of X is the equivalent gain of nonlinearity, 
its value should be larger than zero, as it was proved in Chapter 3. 


For the phase investigating recall the equation (5-19) and re- 





write: 
ja) at ‘=—hb ) mab -- e 
=—=—- z ! Lo te a m fu; @ - i, f (5-23) 
or: 
be a A 
— yo 2 " ec oz $ 4 5 rr 
cos st i = = 4 - . ( o 6h! 
NOG ‘ale Fi -pPi~. oT 
hence: 
O = ein! ——_ : 
wet - . i - ‘ 96 
rf ~ { a a Jj r Ww ee (5 ew) 
in which the value of; 
. f / f J ; E 2: } 
f eae Tan ies WwW a, . ( Ww — /) (5-265) 


Equation (5-25) is shown that the value of @ is a function of 
and the amplitude of first harmonic of error signal, ir the system 


constants is fixed. 
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If we use the equivalent gain for the system investigation, and 


with a result: 


Ge = CTT J A Fig, 0, (5-27) 
i t # es 7G J Le 
Or. 


Q. a AT tat, * 7 od 36, ) 
{ —_ a ’ ) 24 ~ J = 
f - Tg CULT Fac, ¢ Pace fe | (5 28) 


Where om is the amplitude of first harmonic of error signal of forced 
oscillation. 

Equation (5-28) is similar to a linear system, the techniques for 
linear system are applicable. 

5-4, Analog Computer Analysis for a Higher Order System: 

Consider a feedback control system with a transfer function and a 


nonlinearity are shown in Fig. 5-2: 


[= “a | | oO, 


~. 
(y | wa é 7 — ‘ : é r } 
1 ; ti Psy es 
m MX : >| “1 


| | SUS rolySt 2705 Hi 
Las a DOM 
a!) ee i. j 


Fig. 5-2 Block Diagram for a Typical 
Feedback Control System 





in which the restoring force function is: 
E(E) = 2E + 0.020E> (5-29) 
The characteristic curve of restoring force function of nonlinear- 
ity from the calculation is shown in Fig. 5-3. Where: 


E= 1.00; 2.00; 3.00; 4.00; 5.00; 6.00; 7.00; 7.50; 


£(E) 


tt 
— 
e 
oO 
co 


3.84: 5.46: .6:72; 7.50: 7268. .7.15: “6260- 
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For the simulation of the restoring force function of nonlinearity 
from the Computer, it is approximated hy three straight lines, one of 
which with a slope of 2; one of which with a slope of 1; and the other 
with a slope of zero. The circuit setup in the computer is shown in 
Fig. 5-4, all components used are schematically in the diagram, 

The characteristic curve of restoring force function from the com- 
puter is shown in Fig, 5-5; (a) shows the waveforms of the input and the 


output; (h) shows the characteristic curve of the nonlinearity. 


o qT, tb dy b Ks 
Set Pi 
i { 
yi al - fa ! E.; 
A | vp 
| Kw adit 
-— Re Oe | 
a 5) ' ee F me = 
amereeememe fb yy ~ m F : a 
mit = ib 


Fig. 5-4 Analog Computer Setup for 
€(E) = 2E - 0.02083 Simulation 


The operatioral block diagram and the operational circuit diagram 
are shown in Fig. 5-6 and Fig. 5-7. The all components used in the 
operational circuit diagram are sechematically in the diagram, 

The results from the computer are condensed in Table 5-1, and with 
response curves shown in the ¢vwu plane and ¢ 4 plane are shown in Fig. 5-5 


and Fig. 5-9 respectively. 


96 














Je 
| 


. a) i 


1 
| 
: 


‘1 
; Se SS Sy 








| - 
\ 
i 
~ | 
aaa | wy) 
Al) milters!) |. = 
| ' 
> 
| | 
a = i | 
: | | 7 


en eet —— A a —— 
' | ! ' 
, i 
| 
x } 
' } 
=> one = Fe rt —_ = 


ial | = 
a tb 
= -_ a 3 eee oe Mee eek 
~, < | 


oe oe oe eS eS = I 
ge 


= 


=-=<= 
== a a a <= 
4 
q 


J 
i= Coho oe 
| . - 


_——— i Lt FEL 
reba! Bt 
Tas 


rc On 


fn 
7 


OCR 
A FF 
> O 
\ o- Lz 
; 4 Sy 
oa 
; -!> Oo 
oe iS) 
gs rn 
1 ©. 9 Es 
) cs oS 
: o~ 
_| wy Cr 
ed 
wy Ge 
“c- 
. a 
| 
\ 
| 


Se f° 6 EF ES ©” EE EEE = — = 








OGOOOF?P—— 








Operational Block Diagram setup 


Fig. 5-0 
in the Analog Computer 
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Table 5-1. Data for Response Curve from Computer for a Higher Order 
System; in which €(E) » 2E - 0.020E3, G, ({t) = 3.0 Cos@t +8) 


(a) Frequency Increasing: 


Frequency of Amplitude of Phase of Amplitude of Phase of 


Input Error Signal Error Out put Out put 
(i ) ( &) (c) (@) (e) 
0.314 0.40 270.00 3.10 0.00 
0.377 0.55 270.00 3, 30 0.00 
0.440 Oey as) 262.00 3.40 0.006 
0.503 0.90 255.00 3.50 3.80 
0.565 1.10 246.00 3.60 9.00 
0.628 1.25 241.00 3.80 12.50 
0.755 }.70 235.00 4.30 20.00 
0.880 2.50 220.00 5.00 24.50 
1.000 7.40 105.00 6. 30 107.50 
1.130 6.90 64.00 5.10 120.00 ° 
1.260 6.50 37.50 4,25 137.50 
1. 380 6.10 25.00 BU 155.00 
1.510 5.60 15,00 3.00 165.00 
1.640 Sit) 10.00 2550 172.00 
1.760 4.50 7.50 2.00 180.00 
1,890 4.50 5.00 1.70 186.00 
2.200 4.00 0.00 1.20 195.00 
2.500 3.50 0.00 0.85 210.00 
2,810 5.30 0.00 0.60 220.00 
3.140 3.20 0.00 0.40 232.00 
3.460 3.10 0.00 0.30 245.00 
3.770 3.00 0.00 0,25 260.00 
4.090 3.00 0.00 0,20 270.00 


(b) Frequency Decreasing 


Frequency of Amplitude of Phase of Amplitude of Phase of 


100 


Input Error Signal Error Output Outrut 
(a) ( c,) (ce ) Ch) (Cc ) 
2.510 3.50 0.00 0.78 210.00 
2.200 4.00 0.00 1.20 195.00 
1.890 4.50 5.00 1.65 185.00 
1.760 4.90 T50 2.00 180.00 
1.640 Ep As) 10.00 2.50 171.00 
1.510 5.70 15.00 3.00 162.00 
1.380 6.10 25.00 3.50 155.00 
260 6.60 37.50 4.25 137.00 
1.130 6.60 64.06 5.00 121.00 
1,000 7.40 105.00 6.25 106,00 
0.943 / / 6.80 78.00 
0.880 6.70 150.00 6.00 25.00 
0.754 1.80 235.00 4,30 20,00 
0.628 LaZ> 248.00 3.80 16.00 
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CHAPTER VI 
CONC LUSTONS 

The aim of this dissertation was to present a method for investigat- 
ing the conditions for existence of the forced oscillations of a nonlinear 
feedback control system. This method is called "harmonic linearization 
method", it is based on the harmonic balance and iteration method, From 
this method, the conditions for existence of the fundamental frequency 
and the sub-harmonic of order 2 and 3 are investigated, and also a equi- 
valent gain of the nonlinearity in terms of input amplitude and frequency 
has been developed, 

The term of the ''response curve" as used here is defined as the combina- 
tion of two components, one is the attenuation frequency; it is either 
shown in the ''Error vs Frequency" plane or "Output vs Frequency" plane, 
and the other is the phase frequency response, it is with a difference 
value from the different representation of attenuation frequency response. 

The results of this analysis show that a nonlinear feedback control 
system may possess more than one type of forced oscillations, which type 
will exist, it is depending on the characteristic of the system and the 
characteristic of input forcing function. 

For any nonlinear feedback control system, the response curve of 
either fundamental or sub-harmonic forced oscillations, it is usually a 
sudden jump in amplitude and phase when a jump transient take place. The 
range of jump frequency depends on the value of damping of the system and 
the amplitude of the input; usually it is increased as the damping decreas- 
ing, and decreased as the input amplitude decreasing. A vertical tangents 
locus investigating for calculating the jump frequencies in terms of system 


constants and the amplitude of input was investigated. 
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Due to the resonance phenomena and the linear approximations of re- 
storing force function, there is some difference between the theoretical 
and the experimental response curves, It ifs only existant in the jump 
frequency range, other than this range of frequencies, it is very close 
to each other. 

Sub-harmonic oscillations are undesired for a nonlinear feedback 
control system, it is actually an unstable phenomenon. Under suitable 
conditions, a multiple order sub-harmonic forced oscillation may exist 
in a system with a constant input forcing function. There ts also 2 
jump phenomena for each order of sub-harmonic oscillations. The condi- 
tions for the existence of the 2nd and 3rd order sub-harmonic forced 
oscillations are investigated, and also theoretical response curve equa- 
tions for each case are developed. 

It is possible an important advantage that the design and analysis 
of a nonlinear Feedback control system by using the harmonic lLineariza-~ 
tion method, by means for invest igating an equivalent gain of the non- 
linearity. In this case, all principles and techniques applied for linear 


system are applicable. 
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